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Axion Cosmology and the Energy Scale of Inflation 
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We survey observational constraints on the parameter space of inflation and axions and map out 
two allowed windows: the classic window and the inflationary anthropic window. The cosmology of 
the latter is particularly interesting; inflationary axion cosmology predicts the existence of isocurva- 
ture fluctuations in the CMB, with an amplitude that grows with both the energy scale of inflation 
and the fraction of dark matter in axions. Statistical arguments favor a substantial value for the 
latter, and so current bounds on isocurvature fluctuations imply tight constraints on inflation. For 
example, an axion Peccei-Quinn scale of 10^^ GeV excludes any inflation model with energy scale 
> 3.8 X lO'^* GeV (r > 2 x 10~^) at 95% confidence, and so implies negligible gravitational waves 
from inflation, but suggests appreciable isocurvature fluctuations. 



I. INTRODUCTION 

Early universe inflation and the QCD axion provide 
explanations for otherwise mysterious features of the uni- 
verse. Here we argue that assuming both at once leads to 
very significant constraints on their central parameters, 
and to highly falsifiable predictions. 



A. Energy Scale of Inflation 

Inflation is the leading paradigm for early universe phe- 
nomenology ^-[3]. Its mechanism and the values of its 
central parameters are unknown, however. One central 
parameter is the energy scale of inflation Ei, defined as 
the fourth root of the inflationary potential energy den- 
sity, evaluated when the modes that re-enter the horizon 
today left the horizon during infiation. Ei is subject to 
both theoretical and observational constraints, as illus- 
trated in Figs, [l] & [2| 

A multitude of inflation models involving a broad range 
of energy scales have been discussed in the literature, in- 
cluding chaotic inflation [1] , brane inflation O |S] and oth- 
ers [7]. However, very high Ei has been argued to be the- 
oretically problematic, at least for single field slow-roll in- 
flation, because it involves super-Planckian displacements 
of the infiaton field |1] . From Ref . 9J the region in which 
(j) moves at least two Planck masses is £' > 2.4 x 10^^ GeV. 
The intuition that high Ej is problematic seems borne 
out in many string theory models; an example is D-brane 
models [H]. 

Very low Ei has been argued to be theoretically prob- 
lematic also. Naive consideration of families of poten- 
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tial energy functions suggests that Ej < 2 x 10^^ GeV 
(r < 0.01) is non-generic [10]. One of the most strik- 
ing successes of high Ej potentials is that they can natu- 
rally predict ris ^ 0.96, and generic low-energy potentials 
fail to make this prediction. Low Ej potentials have the 
slow-roll parameter e exponentially small, so that the ob- 
servation Us ^ I - 6e + 2t] = 0.960 ± 0.013 [U imphes 
rj — — 0.02±0.0065, leaving us wondering why rj is so small 
when it could just as well have been of order —1. This 
problem is not alleviated by anthropic considerations [12] . 
By using the observed value of density fiuctuations, and 
setting e < 10~^ as the boundary, £■ < 6.7 x 10^^ GeV de- 
fines this problematic low-scale region. These theoretical 
issues for inflation are indicated by the vertical regions 
in Fig. [T] Although there are inflation models in the lit- 
erature at energy scales both above and below this naive 
window, the debate about whether they are generic con- 
tinues. 

With theory in limbo, we turn to observational guid- 
ance. High Ej implies a large amplitude for primor- 
dial gravitational waves (GWs). Ej > 3.8 x 10^^ GeV 
(r > 0.22) is ruled out by WMAP5 plus BAO and SN 
data [TT], as indicated by the orange region of Fig. [2J 
Possible future searches for primordial GWs have rightly 
been a focus of attention. In this article we emphasize the 
additional information that can be learned from isocur- 
vature fiuctuations. 



B. Axion Physics 

The QCD Lagrangian accommodates a gauge invari- 
ant, Lorentz invariant, renormalizable term oc ^E^-B", 
with 6g [~tt,tt], that manifestly breaks P and T symme- 
try. Precision bounds on the electric dipole moment of 
the neutron constrain \9\ < 10~^°. The striking smallness 
of this parameter, which the standard model leaves unex- 
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FIG. 1: Naive expectations for the energy scale of inflation 
Ei and the axion PQ scale /„. For Ei > 2.4 x 10^'' GeV, 
the inflaton must undergo super-Planckian excursions in field 
space (in single field models). For Ei < 6.7x 10^* GeV, generic 
infiation potentials fail to reproduce the observed nearly scale 
invariant power spectrum. For fa > 2.4 x 10^* GeV, the PQ 
breaking scale is super-Planckian. For fa < 10^^ GeV (and 
fa ^TeV), the PQ breaking is in the "desert" of particle 
physics and non-trivial to achieve in string theory. This leaves 
the region labeled "naive window" . 



plained, defines the strong P and T problem (a.k.a. "CP 
problem"). After introducing a new asymptotic (or al- 
ternatively, classical) Peccei-Quinn (PQ) symmetry [13] 
which is spontaneously broken, the effective 6 becomes a 
dynamical variable, and relaxes toward extremely small 
values. The consequent approximate Nambu-Goldstone 
boson is the axion [14, IS^. 

The simplest axion models contain only one phe- 
nomenologically significant parameter: fa, the scale at 
which the PQ symmetry breaks. The zero temperature 
Lagrangian for the complex field cj) = pe'^/\/2 is 



C 



:faid,ef + Jd,pf^2A''smH0/2)-Xm'-fam' 



(A w 78MeV, p is irrelevant at low energies). Acceler- 
ator bounds require fa to be well above the electroweak 
scale, and stellar astrophysics constraints place consider- 
ably higher limits. Given that electroweak values for /„ 
are ruled out, economy suggests that fa could be asso- 
ciated with unification or Planck scales, rather than the 
"desert" of particle physics or super-Planckian scales, as 
indicated by the horizontal regions in Fig. [T] This intu- 
ition for fa seems borne out in string theory, where fa 
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FIG. 2: Observational constraints on the energy scale of infla- 
tion Ei and the axion PQ scale fa are shown in the top (bot- 
tom) panel for inefficient (efficient) thermalization at the end 
of inflation. The thick red diagonal line is fa = Max{TGH ~ 
Hj/2'JT, r„,ax — ectt Ej}, (toff « in top, e^a = 10~^ in bottom, 
with Eotf = 10~^'^, 1 indicated). Above this line is the infla- 
tionary anthropic scenario and below this line is the classic 
scenario. The region in which there is too much isocurvature, 
aa > 0.072, depends on the axion fraction Ra = ^h/Ccdm of 
the CDM; the purple region applies for any Ra, the blue region 
is for Ra > 0.25% (which is expected at 95% confidence), and 
the cyan region is for Ra = 100%. The green region has too 
much axion CDM: ^a > 2.9 eV. Each constraint is divided into 
two parts: the darker part is for a conservative value x = 1/20 
and the lighter part is for a moderate value x = 1. The or- 
ange region has excessive GWs amplitude: Qt > 9.3 x lO"''. 
The yellow region has too much axion interaction in stars 
(darker is firmly ruled out, lighter is for some analyzes). The 
brown region is excluded by the laboratory ADMX search. 
The dashed cyan, orange, and brown lines are future targets 
for isocurvature, GWs, and ADMX searches, respectively. 



typically lies at or just above the GUT scale, and much here: 
lower values are non-trivial to achieve |16j.^ Such high 
values of fa correspond to large contributions from axions 
to cold dark matter (CDM). Indeed, it is only after selec- 
tion effects are taken into account that the ratio of axion aa 
density to entropy is small enough to be consistent with 
observations [HI [19]. When these effects are included, 
one finds that the expected density of dark matter in 
axions is close to the amount of dark matter actually ob- 
served [20j. where 
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C. Cosmological Observables, Summary 



Quantum fluctuations in an effective inflaton field give 
rise to the standard adiabatic fluctuations that have 
grown into our cosmologically observed large-scale struc- 
ture. If the PQ symmetry undergoes spontaneous sym- 
metry breaking before the end of inflation, quantum fluc- 
tuations in the consequent light axion field give rise to 
isocurvature fluctuations. The amplitude of the isocur- 
vature fluctuations grows with S/, so upper bounds on 
the amplitude of isocurvature fluctuations imply upper 
bounds on Ej. 

The purpose of this article is to delineate these bounds, 
extending earlier work on this subject such as [T71 \7T\ - 
[25] . The bounds depend sensitively on the fraction Ra 
of CDM in the form of axions, which in turn depends 
not only on fa , but also on the local initial misalignment 
angle 9i e [— tt, vr]. These constraints are shown in Figs.[2J 
for different choices of the axion CDM fraction. We will 
estimate this fraction using statistical arguments. 

In Section [n] we calculate the production and late time 
abundance S.a of axions and the amplitude aa of isocurva- 
ture fluctuations, as well as reviewing the amplitude Qt of 
primordial GWs. These three observables depend on two 
micro-physical parameters: the PQ scale fa (or equiva- 
lently, the T = axion mass rua) and the energy scale 
of inflation Ej (or equivalently, the Hubble scale of in- 
flation Hj), and on one "environmental" parameter: the 
misalignment angle 9i. We summarize the final formulae 



^ For example, in weakly coupled heterotic string theory, the 
model-independent axion has its PQ scale given by fa = 
aurhpi/{2TT\/2). A unified coupling au = 1/25 then gives 
fa ~ 1.1 X lOifSGeV [TBI [IT]- 
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(definitions are given below). 

The most recent observational bounds from WMAP5 
combined with other data are |11] 

^a<2.9eV, aa< 0.072, Q* < 9.3 x 10-^ (6) 

thereby constraining the two micro-physical parameters 
Ei and fa- 

These expressions for ^a and aa only apply if the PQ 
symmetry undergoes spontaneous symmetry breaking be- 
fore inflation and is not restored thereafter. This is true 
if fa exceeds the Gibbons-Hawking temperature during 
inflation and the maximum post-inflationary thermaliza- 
tion temperature, as we discuss in Section [lI C| Inefficient 
thermalization leads to constraints displayed in Fig. [2] 
(top), while efficient thermalization leads to constraints 
displayed in Fig. [2] (bottom). 

In Section [lll[ we present a statistical estimate for the 
axion abundance to place additional constraints on out 
parameter space. In Section |IV[ we conclude by dis- 
cussing the implications for infiationary model building 
and future prospects. 



II. AXION COSMOLOGY 

In this section, we review the production of axions in 
the early universe, their abundance in the late universe, 
and the amplitude of isocurvature fluctuations, following 
Refs. [T71 [23], and explain and derive eqs. ([T)~([5]). We 
focus on axion production from the so called "vacuum 
misalignment" mechanism only. This provides the most 
conservative constraints. Additional production mecha- 
nisms, such as cosmic string decay, are subject to larger 
theoretical uncertainties (e.g., see [H]). 



A. Onset of Axion Production 

In an expanding flat FRW background at temperature 
T with Hubble parameter H{T), the phase field 9 of bro- 
ken PQ symmetry satisfies the equation of motion 



+ 3H{T)9 



a2 






V{9,T), 



(7) 



where dots indicate derivatives with respect to co- 
ordinate time. Here V{9,T) is the temperature depen- 
dent potential induced by QCD instantons. At zero tem- 
perature, V{9,0) = A'^(l - cos6'), where A « 78MeV 
sets the scale of the vacuum energy of QCD.^ For small 
values of the axion field, the potential is approximately 
harmonic: 
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l/(0,T)«-m,(T)Va' 



(8) 



The mass is temperature dependent, with high and low 
T limits given by 



ma{T) sa ma{0) 



i^r 



for T > A 
for T < Ar 



Qi 



(9) 



where Aq ^ 200 MeV is the scale at which QCD becomes 
strongly coupled, b = 0(10"^) depends on detailed QCD 
physics, and ma(0) is the zero temperature axion mass, 
related to the PQ scale fa and A by niaiO) — h? / fa-^ 
The temperature dependence is illustrated in Fig. |3] 

In the early universe, the axion field is effectively mass- 
less, and so the right hand side of eq. ^ is negligible. 
Hence the zero mode of the axion field is essentially frozen 
due to Hubble friction. When the temperature T drops 
below T„,„, defined by 



3H(T„,J « m,(r_ 



(10) 



the axion field will begin to oscillate, producing axions. 
Since this occurs during the radiation dominated era, we 
have^ 
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FIG. 3: Axion mass m^ (solid) and axion energy density pa 
(dashed), as a function of temperature T (or time t). The 
red and blue curves are for PQ breaking scales fa ~ lO'^^ GeV 
and fa ~ 10^^ GeV, respectively. The arrow indicates the 
effect of increasing fa- The decreasing orange curve is the 
Hubble parameter 3-H". We have taken {6'^) = tt^/S, although 
this is modified by the intervention of inflation, as explained 
in the text. For simplicity, we have here only kept track of 
the variation with temperature in the number of relativistic 
degrees of freedom before (g-t — g,s = 61.75) and after (g* — 
g,s ~ 10.75) the QCD phase transition, and taken x = 1- 



where rhpi w 2.4 x 10^® GeV and g^. (T), the effective num- 
ber of relativistic degrees of freedom, depends on whether 
T„3^ occurs before or after the QCD phase transition: 
5,(r„,J = 61.75 for T„,, > Aq and 5*(r„,J = 10.75 for 
Tasc < Aq. Eqs. (|9 11 1 allow us to solve for T„^„ in terms 
of /a: 



T„ 



0.36 Aq/^ A 



1/3 frnm\ 
[fa J 

0.55A(^) 
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for fa<fa, 
for fa> fa, 



(12) 
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(11) where fa, which reflects the break in eq. ([9|, is defined in 
eq. (16) below. 
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A is set by Aq and quark masses: A^ = 
0.56. 

If there are A^ distinct vacua for d, then we should replace fa 
by fa/N here and throughout the article. However, any Af > 1 
models are expected to have a large overabundance of energy 
density from domain walls, unless inflation intervenes. 
We assume here that the universe before BBN is adequately de- 
scribed by conventional physics. See Ref. 27 for other scenarios. 



B. Density of Axions 

At the onset of production (when T = T^^f) the axion 
energy density is 

Pa{T...) ^\ma{T^..rf!{9')fie^)x. (13) 



Here (6*^) is the spatial average over our Hubble volume 
of the square of the initial misalignment. In terms of its 
mean 9i and standard deviation ag, {9"^) = Of + ag. If the 
axion field is established before (or during the early stages 
of) inflation, then spatial variations in are smoothed out 
over our Hubble volume {V'^O/a? — >• 0). Then 9i = (0) 
in our Hubble volume is an angle drawn from a uniform 
distribution: 0i G [— 7r,7r], with a small variance that we 
discuss in the next subsection. On the other hand, if the 
axion field is established after inflation, then 6i — (0) — 0, 
with variance ag — tt'^/S due to small scale variations. 
f{0l) is a fudge factor acknowledging anharmonicity in 
the axion potential; for 0i — >■ 0, f{0j) — > 1. Finally, x 
is a dimensionless correction factor due to temperature 
dependence during formation. In our numerics, we take 
Aq = 200 MeV, h = 0.018, and absorb all theoretical 
uncertainties into %. A conservative value is x = 1/20 
and a more moderate value is x = 1;^ both values are 
reported in Fig. 2. 

By following the redshift as the universe expands from 
axion-formation to today, we can convert this initial en- 
ergy density into a prediction for the present density, as 
illustrated in Fig. [3] Since we are focusing on the zero- 
mode, the axions form a non-relativistic Bose-Einstein 
condensate. At late times (say today's temperature Tq), 
the axion energy per photon is 



Ca 



PajTo) _ ma(To) Pa(ro3o) sJTq) 



(14) 



where we have exploited the fact that entropy density 
s{T) in a comoving volume is conserved. The entropy 
density is given by 



s{T) 



2^ 



a*s{T)T'' 



(15) 



where g*s{T) is the effective number of relativistic degrees 
of freedom for s{T) |30j. Note that g^siT^..) = 5*(r„,J, 
g*s{To) comes from photons and neutrinos: g*s(To) — 2 + 
I x6x YY = 3.91, and r?T.o(To) « TOa(O). Since the number 
density of photons nj{To) = 2C,{2>)Tq / -k'^ depends on tem- 
perature in the same way that the late time axion energy 
density PaC^o) does, the quantity S,a is a temperature- 
independent, or equivalently time-independent, measure 
of the axion abundance. In contrast, the commonly used 
quantities fi^ and h?^a do not tell us anything fundamen- 
tal about our universe, since like T, they are effectively 
alternative time variables that evolve as our universe ex- 
pands. These different measures of axion density are re- 
lated by Vtah'^ « 0.0019 (^Jl eV)(To/l K)^, which at the 



present epoch (Tq = 2.725 K) reduces to Vlah? « ^a/26eV 
{h is the dimensionless Hubble parameter). 

Combining all this information yields eq. (IT]). Note that 
higher /^ (for a fixed value of {9'^)) corresponds to higher 
axion energy density, as seen in Fig. [3] The reason for 
this is that higher fa corresponds to lower ?Tia(0), so that 
the onset of axion production, when 3H{T) has fallen to 
ma{T), occurs later. Hence there is less redshifting of 
the axion energy density after production (furthermore, 
Pa{T,,,) ishigher if /a < fa). 

We locate the boundary between the low and high fa 
limits by equating the two expressions for ^q . They match 
when fa^ fa, where 



fa = 0.26 



(16) 



The observed total density of cold dark matter in our 
universe ^cdm ~ 2.9 eV implies £,a < 2.9 eV [TT] . 



C. Fluctuations from Inflation 

During inflation, the universe undergoes an approxi- 
mately de Sitter phase with Hubble parameter Hj. Quan- 
tum fluctuations during this phase induce several kinds 
of cosmological fluctuations. 

• Adiabatic density fluctuations are generated, with 
an approximately scale-invariant spectrum. The 
measured amplitude is Q sa 1.98 x 10~^ [II]-^ 

• Primordial gravity waves are generated, with an ap- 
proximately scale-invariant spectrum whose ampli- 
tude is given in eq. ^. WMAP5 plus BAO and 
SN data imply the bound: Qt < 9.3 x 10"^ (95%). 
Since Q is measured and Qt is set by E'/, the ten- 
sor to scalar ratio r = {Qt/Q)^ is often used to 
characterize the scale of inflation. Using the same 
notation as [TT], it is bounded by r < 0.22 (95%). 

• Any other light scalar fields, such as the axion, are 
imprinted with fluctuations during inflation, simi- 
larly to gravitons. The power spectrum of a canon- 
ically normalized scalar field 0, such as (f)a, in de 
Sitter space has a scale-invariant spectrum (e.g., see 

m) 
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(17) 



^ See Refs. 1281 1291 for precise estimates of axion abundance. 



S Here Q = f AB(fc = 0.002 Mpc"!) of Ref. [Tl]. 



It is essentially a thermal spectrum at Gibbons- 
Hawking tcmperatm'e Tgh — Hi/2tt. Fluctua- 
tions in the misalignment angle in fc-space are 
scaled as ag = cTa/fa, per eq. ([5 1. We write the 
corresponding fluctuations in real space as aa == 
7-ff//27r, where 7 = 0(1) is a dimensionless con- 
stant. Ref. [11] effectively takes 7 = 1 and Ref. [22] 
argues that observations are sensitive to length 
scales corresponding to 7 ~ 4, while in our figures 
we have taken a moderate value of 7 — 2. These 
fluctuations provide a lower bound on ^^ and, as 
we discuss in the next subsection, on isocurvature 
fluctuations. 

In the preceding discussion, we assumed the existence 
of a light axion field during inflation. This is true only if 
PQ symmetry is broken before inflation. Furthermore, if 
PQ symmetry is restored after inflation, the fluctuations 
will be washed out. PQ symmetry can be restored either 
by the Gibbons-Hawking temperature during inflation, or 
by the maximum thermalization temperature after infla- 
tion T^^^7 To characterize the maximum thermalization 
temperature, we use a dimensionless efficiency parame- 
ter e„ff defined as T^^^ = e^nEj, with < e^^ < 1, with 
eeff ^ 1 expected. 

A robust criterion for the presence of the axion during 
inflation with fluctuations that survive is 



/a>Max{ToH,T_J. 



(18) 



If this condition is satisfied, inflationary expansion im- 
plies that 6i e [— 7r,7r] is drawn from a uniform distri- 
bution. By postulating that 9i is atypically small in our 
neighborhood (i.e., in our Hubble volume) one can accom- 
modate large fa ■ This defines what we term the anthropic 
regime (see Fig.[2|. 

Alternatively, if fa < MaxjTGH, r„^^}, then either 
there is no axion during inflation or its effects are washed 
out after inflation. In this case 6^ fluctuates throughout 
our observable universe, with variance 7r^/3, and there are 
no appreciable axion-induced isocurvature fiuctuations. 
This defines what we term the classic regime (see Fig. [2]). 



D. Isocurvature Fluctuations 

Fluctuations in the local equation of state 6{ni/s) ^ 
at fixed total energy density Sp = arc known as 



isocurvature fluctuations. (In contrast, fluctuations with 
S{ni/s) = and 6p ^ are known as adiabatic fluctua- 
tions.) Since the axion is essentially massless in the early 
universe, at temperatures much greater than the QCD 
phase transition (T >> Aq), its energy density is entirely 
negligible at early times. Hence, at such early times, fluc- 
tuations in the number density of axions (established by 
de Sitter fluctuations during inflation) do not alter the 
energy density of the universe. Later, for temperatures 
below the QCD phase transition {T < Aq), the axion ac- 
quires a mass and a significant energy density (see Fig.lsl), 
but any such fluctuations cannot alter the total energy 
density of the universe, by local conservation of energy. 
In the early radiation dominated era, this means that 
fiuctuations in the axion energy density are compensated 
by fiuctuations in photons and other relativistic fields. 
Hence, these are isocurvature fluctuations.® 

To quantify the amplitude of isocurvature fluctuations, 
it is useful to introduce the fractional change in the num- 
ber density to entropy density ratio: 



S,^ 



rii/s 



Hi 
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ST 



(19) 



For adiabatic fiuctuations. Si = 0. We assume that this 
is true for all species other than the axion. Isocurvature 
fluctuations in the total energy density involve a sum over 
all massive species and radiation: 



= (5/3i3o = TTlaSUa + /. miSrii 



Apr 



5T 
Y' 



(20) 



These two equations will be used to obtain an expression 
for the corresponding temperature fluctuations. 

Initially the energy density of the axion field is a small 
fraction of the ambient total density, so eq. (20) gives 
5T /T ^ 5na/na, and Sa — Sna/ua- Since Ua oc 9^ (ig- 
noring anharmonic effects), this implies 



Sa 



(21) 



Assuming 59 = 6 — {6) is Gaussian distributed^, we can 
calculate (5^) in terms of 9i = {6) and a^ = {{69)'^), as 



{sD = 



2al{2e^,+al) 



r2A2 



(22) 



^ The maximum thermalization temperature should not be con- 
fused with the reheating temperature, which can be somewhat 
lower [32| . The maximum thermalization temperature is the 
maximum temperature of the thermal bath post-inflation, while 
the reheating temperature is the temperature at the end of the 
reheating phase, i.e., at the beginning of the radiation era. 



° Later, around the onset of the matter dominated era, these 
isocurvature fiuctuations are converted to adiabatic fluctuations, 
responsible for the familiar gravitational structures in our uni- 
verse. 

^ This is a good assumption in the regime where the axions com- 
prise a significant fraction of the dark matter, i.e., 9^ 2> o-g- 
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Note that if Of < ct^ then {SD = 2, while if 9f > ct^ 

then(52^=4a2/^f. 

The most important axion induced temperature fluctu- 
ations are those on the largest scales. Such fluctuations 
enter the horizon well into the matter dominated era, 
where pr can be ignored. This implies^'^ 



6T\ 






"-'a: 



(23) 



where ^m is the total matter energy density per photon, 
whose measured value is ^„i = 3.5 eV [llj. 

Following [TTj , we define aa to be the fractional contri- 
bution to the CMB temperature power spectrum due to 
axion isocurvature: 



_ {{5T/T)l 



iiST/T) 



2 ) 

tot/ 



(24) 



Using the relationship between Sa and ST/T in eq. ( 23 ) 
and the preceding expression for (S*^) we obtain 



25 {{ST/T)L) {9^ + aly ' 



(25) 



COBE measured (and WMAP confirmed) the root-mean- 
square total temperature fluctuations to be {{ST/T)'^^^) sa 
(1.1 X 10"^)^, averaged over the first few I. Using the 
expression for S,a given in eq. jTl) we can write the isocur- 
vature fluctuations as in eq. ( ST. This must be consistent 
with the latest observational bound Ua < 0.072. Here 
we have used ag of Ref. [11] , which assumes isocurvature 
fluctuations are uncorrelated from curvature ones. 



III. DIRECT AND STATISTICAL 
CONSTRAINTS 

It is conceivable that the axion abundance is negligible 
(but see the following subsection). This scenario (case (i)) 
requires 9^ ^ ag. By demanding aa < 0.072 (the current 
isocurvature bound) and using eq. (pLwe obtain the most 
conservative bound, studied in Ref. [TTj, corresponding to 
the purple region marked "Any i?^" in Fig. [2] 

At the other extreme, if axions are the dominant form 
of dark matter in the universe (case ni), then 9f >• ag. 
Again demanding Ua < 0.072 in eq. ^, with 9i deter- 
mined from eq. (fTj) with ^q = ^cdm, the excluded region 
expands to include the cyan region marked "i?a = 100%" 
in Fig. [2] (as well as the blue region marked ^^Ra > 
0.25%"). 



^^ Due to the Sachs-Wolfe effect, there is a 20% enhancement to 
l|23|l, but we will not go into those details here. 



Each of these three regions are bifurcated by a line. In 
all three cases, the rightmost part gives the most robust 
constraint, coming from a conservative value x — 1/20, 
while the leftmost part extend the constraints using a 
more moderate (and more speculative) value % = 1. This 
comes from our uncertainty in the total axion abundance. 



A. Statistics of a Two-Component Model 

The viability of large fa axion cosmology depends on 
taking selection effects seriously, since they produce a 
higher dark matter density ^ than observed in most Hub- 
ble volumes. In particular, the density of a typical galaxy 
scales as p ~ f"*. Taking into account that denser galaxies 
have fewer stable solar systems due to close encounters 
with other stars, etc., it has been found that typical stable 
solar systems in large fa axion models reside in Hubble 
volumes where ^ is comparable to the observed value |20| . 

Here we draw out a statistical implication for the pre- 
dicted axion abundance, if there is a second contributor to 
the dark matter density. Consider the hypothesis that the 
total CDM (^cdm) is comprised of axions (^a) and some 
other component, say WIMPs {£,w)'- Ccdm = £,a+S,w- The 
unknown separate axion and WIMP abundances should 
be drawn from prior distributions determined by under- 
lying micro-physical theories. For axions in the large fa 
regime, above any inflation temperatures, this scenario 
implies that the initial misalignment angle 9i is uniformly 
distributed. In the regime where the axion abundance 
is non-negligible (6*1 ^ a^), but still sufficiently small 
that we can ignore anharmonic effects (6*1 <C 1), we have 
^a oc 6^. Since 9i is uniformly distributed, it is simple to 
show that 



pi±{L) (X 



1 



(26) 



In contrast, we do not have a reliable prior distribution 
P^^li^w) for the WIMP. 

As discussed in Ref. ^\, selection effects depend only 
on the sum ^cdm = ^a + ^w, so the total joint distribution 
for axions and WIMPs is 

As demonstrated in Ref. |20j , the observed value of CDM 
CcDM ~ 2.9 eV is nicely consistent with this distribu- 
tion. Given this, we can focus on the remaining one- 
dimensional distribution for the axion: 



p(ea)«pll(ea)pS(^COM-ea). 



(28) 



IW) 



Unless Pprior is sharply peaked at ^cdmj the axion prior 
(when integrated) disfavors very small values of ^a- For 



example, let us take the WIMP prior to be uniform. We 
can then make a prediction for the axion abundance at, 
say, 95% confidence. Defining ^a implicitly through 



Pda) d^a = 0.05 



(29) 



and solving eq. (29) using eq. (26), we find ^a — 
(0.05)^^cDM = 0.25%^cDM- This says that it statistically 
unlikely - at the 95% level - for axions to comprise less 
than 0.25% of the CDM of the universe (case(ii)). 

By setting ^^ = 0.25% ^cdm, we rule out the blue re- 
gion marked Ra > 0.25% in Fig. [5] with high confidence. 
In other words, without assuming that axions comprise 
all the CDM, we find that on statistical grounds axions 
must comprise at least a non-negligible fraction of the 
universe's CDM, allowing us to extend the excluded re- 
gion in Fig. [2] further towards the upper left. 



B. Additional Constraints 

The preceding applies in the fa > Max{TcH,T,„^^} 
regime, where the initial misalignment angle 0i takes 
on a single constant value in our Hubble volume. For 
fa < Max{TGH,TL,ax}, the misalignment angle varies on 



cosmologically small scales, with average 



73. In 



this regime the isocurvature fluctuations are negligible. 
In this case, bounds arise from the requirement that the 
axion abundance is not greater than the observed total 
CDM abundance: £_a < £cdm ~ 2.9 eV. Using the upper 
expression for S,a in eq. ([l]), with 9f + ag ^- {6"^) = 7r^/3, 
we find that /„ > 2.3 x 10" x^^'"^ GeV is ruled out. For 
the conservative value x = V^O; this excludes the up- 
per part of the green region marked "Too much CDM" 
in Fig. [2j and for the moderate value x — ^ this extends 
the exclusion to the lower part of the green region. ^^ 

Also, ma(0) > 10'^ /zeV is firmly ruled out (and 
m£i(0) > 10* /ieV for some analyzes), since in this regime 
the coupling of axions to matter is too large, affecting the 
physics of stars, such as the cooling of red giants and the 
neutrino flux from SN 1987A 33 (yellow region at bot- 
tom of Fig. [2J) Furthermore, the ADMX search for axion 
dark matter in a microwave cavity detector has ruled out 
axions comprising the bulk of the halo dark matter in 
the following mass window: 1.9/xeV< TOa(O) < 3.3^eV 
(brown band in Fig. [2]) for so-called KSVZ axions, and the 
sub-window 1.98/ieV< 771^(0) < 2.17/ieV for so-called 



^^ If fa > MaxjTGH, r,„ax}, there is another region ruled out with 
too much CDM (green region above thick red line in Fig. pi) 



DFSV axions [3U |3S] . The remaining white region is the 
allowed "classic window" . 

In Fig. ^ (top), corresponding to inefficient thermaliza- 
tion (e^ff ~ 0), the boundary between the anthropic and 
classic regimes is fa = Tgh — Hi /2tt. In Fig. [2] (bottom), 
corresponding to efficient thermalization (e„ff = 10~^, 
with e„ff = lO""'^'^, 1 indicated), the boundary between 
the two regimes is fa — T^^^ = e^aEj. Efficient thermal- 
ization thus opens up a larger "classic window" , but the 
"anthropic window" is essentially unaltered. 



C. Effect of Falling Density During Inflation 

In our analysis, we have treated inflation as occurring 
at one rather well-defined Hubble scale. Although this 
is a good approximation in some inflation models, there 
are others giving an appreciable change in H between its 
value (say Hj) when the modes that are now re-entering 
our horizon left the horizon (55 or so e-foldings before 
the end of inflation), and its value (say i?cnd) at the end 
of inflation. This is particularly relevant to Fig. l2] (top), 
since it implies that the boundary between "anthropic" 
and "classic" regimes is blurred, since Tgh = H/2tt is 
evolving. For high scale inflation, H typically changes by 
an amount of order the number of e-foldings, i.e., 0(10^), 
while for low scale inflation models, H typically changes 
very little. 

If we consider iJ„„d ^ Hj, then the PQ symmetry can 
break during inflation. The resulting cosmology could be 
quite interesting with axion dark matter varying appre- 
ciably from one point in our Hubble volume to another, 
but is ruled out since Q ^ 10^'"'. If PQ breaking occurs 
very close to 55 or so e-foldings before the end of infla- 
tion, then 9i can be smoothed out on today's cosmological 
scales and make the analysis "anthropic" . Otherwise, we 
expect the "classic" analysis to apply as usual, providing 
a ruled out green region in Fig. [2J Hence, we expect such 
corrections to the constraints to be reasonably minimal. 



IV. DISCUSSION 

We have surveyed observational constraints on the pa- 
rameter space {Ej,fa} of inflation and axions, finding 
that most of it is excluded, leaving only two allowed re- 
gions that we term classic and anthropic windows. Part of 
the classic window fa ~ 10^^ - 10^^ GeV will be intensely 
explored by the ongoing ADMX experiment. The region 
indicated by the arrow to the horizontal brown lines in 
Fig. [2] to ma(0) = 10 ^eV is expected to be explored by 
the end of ADMX Phase II, and onwards to ma{Q) = 
100 /ieV some years thereafter |36| . In this window, com- 
paratively little can be concluded about the scale of in- 



9 



flation. From Fig. ^ (top), taking ma{0) = 100 fieV and 
assuming ^a > 0.25% ^cdm, we rule out 8.4 x 10^^ GeV 
<Ei< 1.3 X lO^^GeV.i^ From Fig. [2] (bottom), the up- 
per end of this ruled out region is reduced due to efficient 
post-inflation thermalization. Although we can rule out 
a range of low scale inflation models, these conclusions 
are not exceedingly strong. 

On the other hand, a large fa axion has strong impli- 
cations for inflation. According to both Figs. [2] (top) & 
(bottom), if fa = lOi^GeV then Ei > 5.5 x 10^'^ GeV 
(r > 9 X 10~^) is ruled out at 95% confidence for the 
conservative value x = 1/20, and Ei > 2.6 x 10^'* GeV 
(r > 4 X 10"^") is ruled out for the moderate value 
X — I- The geometric mean is Ej > 3.8 x 10^^ GeV 
(r > 2 X 10~^), which is reported in the abstract. This 
is incompatible with many models of inflation, including 
"classic" models with a single slow-rolling scalar field in 
a generic potential. For example, monomial potentials 
V on (j)P predict r — 4p/A'e, where N^. is the number of 
e-foldings of inflation from when it generated our horizon 
scale fluctuations to when it ended. For such models, A^g 
around 50 or 60 is expected, so any reasonable p is ruled 
out, including 0^ chaotic inflation [3] and the stringy N- 
flation 37 and Monodromy [35] models. The same is 
true for exponential potentials V oc exp(— -v/2p0/mpi), 
which predict r = 16p. 

Evidently there is considerable tension between the 
theoretically appealing large fa and high-scale inflation 
scenario (see Fig. [I]) and the observational constraints 
(see Fig. [2]) . Low-scale inflation may be emerging as fa- 
vored from recent work in string theory. If we consider 
the small subspace (see [5UH1T] ') of presently constructed 
string models that both inflate and agree with the ob- 
served values of Q and n^, we are left with models that 
tend to be at rather small energies, typically r < 10~^ for 
D-brane models and various other scenarios such as mod- 
ular inflation |42] . There are also arguments for very low 
r in the simple KKLT framework discussed in Ref. [33]. 
This allows fa ^ 10^® GeV to be marginally consistent 
with present isocurvature bounds. Although it is highly 
premature to conclude that very low energy scale is a 
generic feature of string realizations of inflation, it is in- 
triguing that many string constructions have this feature. 
(See [37l|38] for interesting exceptions.) 

The Planck satellite, CMBPol, and upcoming subor- 

^■^ The quoted lower end of the ruled out region is the geometric 
mean of the conservative and moderate ^a-scenarios. 



bital CMB experiments should probe well beyond the 
current bound on GWs of r < 0.22, perhaps reaching 
r ^ 0.01. This is indicated by an arrow toward the verti- 
cal dashed red line in Fig. [2] If gravity waves are observed 
in this regime, then the PQ scale fa must be in the classic 
window. 

Our considerations emphasize the fundamental impor- 
tance of improving bounds on isocurvature fluctuations. 
For example, an order of magnitude improvement to 
a ^ 0.007 would push the isocurvature bounds to the 
diagonal dashed cyan line in Fig. I2] (We have indicated 
the improvement for the case where axions comprise all 
the CDM: ^^ = Ccdm-) Detection of isocurvature fluctu- 
ations in this regime has three important implications: 

1. It could be interpreted as evidence for the existence 
of the axion fleld, and assuming this: 

2. It would probe low inflation scales Ei far beyond 
the scope of any foreseeable GW measurements. 

3. It would be evidence that we live in a highly atyp- 
ical Hubble volume, i.e., {_E/,/a} must be in the 
anthropic window. 

Isocurvature modes and tensor modes thus provide com- 
plementary constraints on fundamental physics, making 
it fruitful to study dark matter and inflation in a unified 
way. 
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